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Abstract 

G\ ' 

] The rotation of Mercury is presently captured in a 3/2 spin-orbit resonance with the orbital mean motion. The capture mechanism 
is well understood as the re sult of tidal interactions with the Sun combined with planetary perturbations ( Goldreich and Peale, I96d; 
^^^^ ICorreia and Laskai , |2004|) . However, it is now almost certain that Mercury has a liquid core dMargot et all (2007h which should 



induce a contribution of viscous friction at the core-mantle boundary to the spin evolution. According to Peale and BossI d 19771) 
this last effect greatly increases the chances of capture in all spin-orbit resonances, being 100% for the 2/1 resonance, and thus 
preventing the planet from evolving to the presently observed configuration. Here we show that for a given resonance, as the 
chaotic evolution of Mercury's orbit can drive its eccentricity to very low values during the planet's history, any previous capture 
can be destabilized whenever the eccentricity becomes lower than a critical value. In our numerical integrations of 1000 orbits 
of Mercury over 4Gyr, the spin ends 99.8% of the time captured in a spin-orbit resonance, in particular in one of the following 
three configurations: 5/2 (22%), 2/1 (32%) and 3/2 (26%). Although the present 3/2 spin-orbit resonance is not the most probable 
outcome, we also show that the capture probability in this resonance can be increased up to 55% or 73%, if the eccentricity of 
Mercury in the past has descended below the critical values 0.025 or 0.005, respectively. 
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Mercury's present rotation is locked in a 3/2 spin-orbit res- 
jonance , with its spin axis nearly perpendicular to the orbital 
plane (iPettengill and DvceL [l965). The stability of this ro- 
tation results from the solar torque on Mercury's quadrupo- 
lar moment of inertia, combined with an eccentric orbit: the 
axis of minimum moment of inertia is always aligned with 
the direction to the Sun when Mercury is at the perihe- 
lion of its orb it ("Colombo", '1965^ Goldreich and Peale , 1966 ; 
ICounselman and Shapiro, 1970). The initial rotation of Mer- 
pury was presumably faster than today, but tidal dissipation 
along with core-mantle friction brought the planet rotation to 
the present configuration, where capture can occur. However, 
the exact mechanism on how this state initially arose is not com- 
pletely understood 



In their seminal work, iGoldreich and Peale (Il966h have 



shown that since the tidal strength depends on the planet's ro- 
tation rate, it creates an asymmetry in the tidal potential that 
allows capture into spin-orbit resonances. They also computed 
the capture probability into these resonances for a single cross- 
ing, and found that for the present eccentricity value of Mer- 
cury (e - 0.206), and unless one uses an unrealistic tidal model 
with constant torques (which cannot account for the observed 
damping of the planet's libration), the probability of capture 



into the present 3/2 spin-orbit resonance is on the low side, at 
m ost about 7%, whi ch remained somewhat unsatisfactory. 



Goldreich and P eale ( 1967) nevertheless pointed out that the 



probability of capture could be greatly enhanced if a planet has 
a molten core. In 1974, the discov ery of an in trinsic magnetic 
field by the Mariner 10 spacecraft dNess et al lll974), seemed 
to imply the existence of a conducting liquid core and conse- 
quently an increment in the c apture probability in the 3/2 reso- 
nance. However, according to IGoldreich and Peald (Il967h . this 
also inc reases the capture pro bability in all the previous reso- 
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nances. iPeale and BossI (11977 ) indeed remarked that only very 
specific values of the core viscosity allow to avoid the 2/1 reso- 
nance and permit the capture in the 3/2 confi guration. 

More recently. ICorreia and Laskar (l2004 (hereafter denotec 
by Paper I) have shown that as the orbital eccentricity of Mer- 
cury is varying chaotically, from near zero to more than 0.45, 
the capture probability is substantially increased. Indeed, when 
the large eccentricity variation is factored into the capture, the 
rotation rate of the planet can be accelerated, and the 3/2 reso- 
nance could have been crossed many times in the past. Perform- 
ing a statistical study of the past evolutions of Mercury's orbit, 
over 1000 cases, it was demonstrated that capture into the 3/2 
spin-orbit resonant state is in fact, and without the need of spe- 
cific core-mantle effect, the most probable final outcome of the 
planet's evolution, occurring about 55.4% of the time. In con- 
trast, because the eccentricity can decrease to near zero, all res- 
onances except the 1/1 become unstable, allowing the planet to 
escape from resonance. This mechanism suggests that in pres- 
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ence of core-mantle friction the planet can escape to a previous 
capture in the 2/1 or higher order spin-orbit resonances. 

The present paper continues the work started in Paper I. In 
addition to the effects of tides and planetary perturbations, we 
will consider here also the core-mantle friction effect as de- 



scribed by iGoldreich and Peald (Il967h . since the presence of 
a liquid core inside Mercury is now confirmed by radar obser- 
vations ( Margot etal.^ 2007). In the next section we give the 
averaged equations of motion in a suitable form for simula- 
tions of the long-term variations of Mercury's spin, including 
the resonant motion, viscous tidal effects, core-mantle coupling 
and planetary perturbations. In section 3 we discuss the conse- 
quences of each effect into the spin evolution and evaluate the 
capture probabilities in resonance. Finally, in last section we 
perform some numerical simulations to illustrate the different 
effects described in section 3. 

2. Equations of motion 

We will adopt here a model for Mercu ry which is an exten- 
sion of the model from Poincare (fl910l) of a perfect incom- 
pressible and homogeneous liquid core with moments of inertia 
Ac = Be < Cc inside an homogeneous rigid body with moments 
of inertia A„, < B,„ < C,„, supported by the reference frame 
(i,j,k), fixed with respect to the planet's figure. Tidal dissi- 
pation and core-mantle friction dr ive the obliquity close to zero 



(lYodenll997l : ICorreia et alil2003h . Since we are only interested 



here in the study of the final stages of evolution (where capture 
in spin-orbit resonance may occur), we will neglect the effect 
of the small obUquity variations on the equations of motion. 
Moreover, since for a long-term study we are not interested 
in diurnal nutations, we can average over fast rotation angles 
and merge the axis of pr incipal inertia and the axis of rotation 
( Boue and Laskar . 2006h . Therefore, the mantle rotation rate is 
simply given hy o) - + (p, where is the rotation angle and (p 
the precession angle. 

2.1. Precession torque 

The gravitational potenti al 'V generated at a generic poi nt of 
the space r is given by (e. g. iTisserandL 1 1 89 It ISmarn. 1 1 953h : 
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where terms in (R/r)^ were neglected. G is the gravitational 
constant, m the mass of Mercury, Ur - r/r and P2{x) = 
(3x^ - l)/2 is the Legendre polynomial of degree two. When 
interacting with the Sun's mass, m©, the spin of Mercury will 
undergo important changes. The middle term in the above po- 
tential will be responsible for a libration in the spin, while the 
last term causes the spin axis k to precess around K, the nor- 
mal to the orbit. Since we are only interested in the study of 
the long term motion, we will average the potential over the 
rotation angle and the mean anomaly M, after expanding the 



true anomaly v in series of the eccentricity e and mean anomaly. 
However, when the rotation rate w and the mean motion n - M 
are close to resonance {oj ^ pn, for a semi-intege^^ value p), 
we must retain the terms with argument {20 - 2pM) in the ex- 
pansion 



cos(20 - 2v) 



— f^(P^ ^) cos(26l - 2pM) , 



(2) 



p=- 



where a is the semi-major axis of the planet's orbit and the func- 
tion H{p, e) is a power series in e (Tab. [TJ. The exact averaged 
contributions to the spin a re given in a s uitable form for our 
study by expression (15) in Correia (l2006lfl For zero obliquity 
we have: 



dcL) 
dt 



H(p, e) sin 2{0 + ip - pM - vy) , 



(3) 



where vj is the long itude of the perihehon, Cm - Cm/C — 0.55 
(iMargot etal.[ 120071) . and 



3GmQB-A 3 ^B-A 

P = ^ -n- 

^ ' C 2 



2a' 



C 



(4) 



The Mariner 10 flyby of Mercury provided information on 
the internal structure of the planet, though subject to some 
uncertainty. For the gravity field it has been measured J2 - 
(6.0 + 2.0) X 10'^ andC22 = (1.0 ±0.5) x 10~^ jAnderson et al 



19871) . Modeling the interior structure of Mercury, it has been 
estimated for the s tructure constant ^ = C/(mR^) ^ 0.3359 
(ISpohn et alil200lh . Thus, for the moments of inertia we com- 
pute (B-A)/C = 4C22/f ^ 1.2 X 10""*. 

2.2. Tidal torques 

Tidal effects arise from differential and inelastic deforma- 
tions of Mercury due to the gravitational effect of the Sun. Their 
contributions to the spin variations are based on a very gen- 
eral formulatio n of the tidal potential, initiated by George H. 
Darwinl(ll880l) . The distortion of the planet gives rise to a tidal 
potential, 

^Hr, f) = -h^ (^)' (^)' P2 («,. ■ a,) . (5) 

where R is the average radius of the planet and r' the radial 
distance from the planet's center to the Sun. In general, imper- 
fect elasticity wi ll cause the p hase angle of 'V^ to lag behind 



rect elasticity wi ll cause the p hase angle or 'K* to lag behind 
the perturbation ( Kaulal 1964 1). because there is is a time delay 
At between the perturbation of the Sun and the maximal defor- 
mation of Mercury. During that time, the planet rotates by an 
angle wAf, while the Sun also changes its position. Assuming a 
constant time delay allows u s to linearize the tidal potential an d 
simphfy the tidal equations (lMignard[l97lll980l: [Hut [ll98lh . 



' We have retained the use of semi-integers for better comparison with pre- 
vious results. 



^There is a misprint in the sign of <p in lCorrei3 j200d) . 
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Table 1: Coefficients of H(p, e) to e'. The exact expression of these coefficients is given by H(p, e) = i (2) exp(i 2v) exp(i 2pM) dM. 
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For zero obliquity, the averaged contributions to the spin are 
then given by: 
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and - n/S.t. This tidal model is particularly adapted to de- 
scribe the planets behavior in slow rotating regimes (10 ~ n), 
which is the case of Mercury during the spin-orbit resonance 
crossing. In the present work we will adopt ^2 = 0.4 and 
Q - 50, which yields K = 2.2 x lO^'^yr"^. This choice is 
somewhat arbitrary, but based on the parameter values of the 
other terrestrial planets (Goldreich and Soter. 1966) . 

2.3. Core-mantle friction effect 

The Mariner 10 flyby of Mercury revealed the presence of an 
intrinsic magnetic field, which is most likely due to motions in a 
conducting fluid core (for a review see lNesslll978 ). Subsequent 
observations made with Earth-based radar provided strong evi- 
dence that the mantle of M ercury is decoupled from a core that 
is at least partially molten (iMargot et al. , 2007). 



If there is slippage between the liquid core and the man- 
tle, a second source of dissipation of rotational energy results 
from friction occurring at the core-mantle boundary. Indeed, 
because of their different shapes and densities, the core and 
the mantle do not have the same dynamical ellipticity and th e 
two parts tend to precess at different rates ( Poincare . 1910l) . 
This tendency is more or less counteracted by different inter- 
actions produced at their interface: the torque of non-radial in- 
ertial pressure forces of the mantle over the core provoked by 
the non-spherical shape of the interface; the torque of the vis- 
cous friction (or turbulent) between the core and the mantle; 
and the torque of the electromagnetic friction, caused by the in- 
teraction between electrical currents of the core and the bottom 
of the magnetized mantle. The two types of friction torques 
(viscous and electromagnetic) depend on the differential rota- 
tion between the core and the mantle and can be expressed by 
a single effective friction to rque, T dMathews and Guoi 12005^ 
Deleplace and Cardhl 2006): 



T - Crxd : 



(10) 



where k is an effective coupling parameter, Wc the core's ro- 
tation rate and q - CdC - \ — Cm = 0.45 (IMargot et al. , 
According to Mathews and Guol (2005) we may write 
where Rc is the core radius and v is the 
Adopting RclR = 



K ^ 2.62Vvw/i?, 
effective kin ematic viscosity of the core 



0.77 ± 0.04 dSpohnet al.ll200lb and v = 10 ''m^s we com- 



pute for the 3/2 resonance: k - 5 x 10 yr \ T he uncertainty 
over V is very large, according to Lumb and Aldr idge ( 1 99ll) it 
can cover about 13 or ders of magnitude. As in former studies 
on planetary evolution (lYodeiill997tlCorreia and Laskaiil200ll 
2003), we used the sarne value as th e best estimated so far for 
the Earth dOansl 1 1 9721; [Poirieii 1 1 988h . 



3 



Another important consequence of core-mantle friction is to 
tilt the equator of the planet to the o rbital plane, which results 
in a se cular decrease of the obliquity jRochestei , 1976t Correia , 
reinforcing the previous assumption that the obliquity 
remains close to zero during the last evolutionary stages. 

Since the core and the mantle are decoupled, we need a dif- 
ferential equation for each rotation rate, the coupling equation 
being given by the friction torque ( Peale. 2005; .Correia. .2006,) : 



doj CcK dojc 

— — <5 and 

dt Cm dt 



2.4. Planetary perturbations 



k5 . 



(11) 



When considering the perturbations of the other planets, the 
eccentricity of Mercury is no lon ger constant, but u ndergoes 
strong chaotic variations in time dLaskar , Il990[ 11994, 2008). 
These variations can be modeled using the averaging of the 
equations for the motion of the Solar System, that have been 
compared t o recent numerical integrations, with very good 



agreement ( Laskar et al. . 2004allbl) . These equations are ob- 



tained by averaging the equations of motion over the fast angles 
that are the mean longitudes of the planets. The averaging of the 
equation of motion is obtained by expanding the perturbations 
of the Keplerian orbits in Fourier series of the angles, where the 
coefficients themselves are expanded in series of the eccentric- 
ities and inclinations. This averaging process was conducted 
in a very extensive way, up to second order with respect to the 
masses, and through degree five in eccentricity and inclination, 
leading to truncated secular equations of the Solar System of 
the form 



dw 

— — 1 [Fw + 03(w, w) + 05(w, w)] 
dt 



(12) 



where w is the column vector {z\, ■.■,Z'i,(\, .■■,(^), with Zj - 
exp(itiTy), i^j = sin(/j/2) exp(iQy) and vj and Q, respectively 
the longitude of the perihelion and node. The 16x16 matrix 
F is the linear Lagrange-Laplace system, while 03 (w, w) and 
<I)5(w, M>) gather terms of degree 3 and 5 respectively. 

The system of equations thus obtained contains some 150000 
terms, but its main frequencies are now the precession frequen- 
cies of the orbits of the planets. The full system can thus be 
numerically integrated with a very large step-size of 250 years. 
Contributions due to the secular perturbation of the Moon and 
general r elativity are also included (for more detai ls and refer- 
ences see iLaskaii [l990[ 1 1 9961: iLaskar et al.i l2004ah . 

This secular system is then simplified and reduced to about 
50000 terms, after neglecting terms of very small value (ILaskar , 
[1994). Finally, a small correction of the terms appearing in F 
(EqfT2ll. after diagonalization, is performed in order to adjust 
the hn ear frequencies, in a similar way as it was done bv lLaskai 
(11990) . The secular solutio ns are very c l ose to the direct nu- 
merical integration La2004 (^ Laskar et al. i l2004allbl) over about 
35 Myr, the time over which the direct numerical solution itself 
is valid because of the imperfections of the model. It is thus 
legitimate to investigate the diffusion of the orbital motion of 
Mercury over long times using the secular equations. 



Over several millions years, the eccentricity of Mercury 
presents significant variations (FigUli that need to be taken into 
account in the resonance capture process (Paper I). 



3. Spin evolution and capture probabilities 

The evolution of the spin for zero obUquity is completely 
described when we put together the contributions from the dif- 
ferent effects presented in section |2l 



io^P + T - r/c„ 
(he = r/Cc , 



(13) 



where P is the precession torques (EqO), T the tidal torque 
(Eq|6]l and T the friction torque (EqfTOli. In order to better study 
the capture probabilities in spin-orbit resonances, we will adopt 
a change of variables that is vaUd around each resonance p: 



y -6 + <fi- pM - m , 
and thus, 

y — oj — pn — vj and y - 6j — m . 



(14) 



(15) 



In absence of planetary perturbations, for a Keplerian orbit, 
m - Q and thus y - Cj. Instead of the core rotation cjc we 
will also adopt the diff'erential rotation (5 = w - as variable. 
The above system of equations (fT3] l becomes then: 



y = P(y) + T(y) - CcK,„6 - m , 
6 = P{y) + T(y) - k„,6 , 



(16) 



where k,„ - k/c„,. In general, we like to express y as the sum of 
the precession torque P{y) and a global dissipative torque D{y) 
such that 



y = -/?,„ sin 2y + D(y) , 
where 

y6 



A„ = —H(p,e) 

Cm 



and 



D{y) = -Do 



(17) 



(18) 



(19) 



Do, V and /i being "constant" quantities. Indeed, under this 
linearized form we are able to estimate t he capture probabih- 
ties us ing a simple expression derived by iGoldreich and Peale 
(Il966l) : 



P 

^ cap 



1 + 



n n y 
2 Aw ju 



(20) 



where Aw is the maximal amplitude of libration in resonance, 
i.e.. Aw - v/2yS„,. 



3.1. Effect of tides 

Let us consider first the simplified case where the spin of 
the planet is only su bject to the precession and tidal torques 



( Goldreich and Peal3 . fl966 ). Thus, we may use expression (fTTI) 
to express y, where D(f) is given by expression ©I 



D{y) = r = -— 



p - E(e) + - 
n 



(21) 



with E{e) - N(e)/D.(e). The limit solution of the rotation for 
a constant value of the eccentricity is obtained when D{y) = 
0, that is, for y/n = E(e) - p <^ oj/n - E{e). If the planet 
encounters a spin-orbit resonance in the way to the equilibrium 
position, capture may occur with a probabihty computed from 
expression ( l20l i: 



P -1 

-* cap — ^ 



n n 



(22) 



Using the present value of the eccentricity of Mercury and 
Cm = 1 we compute for the p - 3/2 resonance a probability 
of capture of about 7%, which is unsatisfactory given that this 
is the presently observed resonant state. 

3.2. Effect of core-mantle friction 

We now add the effect of core-mantle friction to the effects 
akeady considered in the previous section. In this case we must 
take into account the rotation of the core, and the complete ro- 
tation rate evolution is given by system ( fTSI l with it = 0. 

The general solution for the differential rotation in system 
( fT6b is given by: 



6(t) = So e" 



+ e 



J to 



+ T(y)] dt' 



(23) 



where 5q is the initial value of 6{t) for t - to- In section 12.31 
we estimate for the present rotation of Mercury 1/k„, ^10^ yr, 
which can be seen as the time scale needed to damp the initial 
differential rotation to zero. Thus, for f » 1 //f„, the evolution 
of the differential rotation is given only by the last term is ex- 
pression ( l23T l. It also means that for time intervals Af <k 1 //c„,, 
we can write for t < to + At: 



6(t) ^6o + 



f ydt' . 

J to 



that is. 



6(t) 5o + r - to = r - tco , 



(24) 



(25) 



where = yo- Sq - - pn, and 0)^0 is the value of the core 
rotation for t - to- This approximation is very useful because 
we can express y by means of expression ( [TtT i with 



Diy) = -—Q(e) 



p - E(e) + - 
n 



CcK 

(y - y^J 

Chi 



K 



-Q(e) 
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n n 
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where - CcKn/[KQ.{e)]. It allows us easily to compute the 
capture probabilities using expression ( l20l l provided that we are 
able to estimate correctly y^^ just before the planet crosses the 
resonance: 



cap 



p-E(e)-x 



1 + 



2 Aoj 



-X 



(27) 



According to expression (fTTT i we have ioc - k6. Since k <& 
Aoj, the core is unable to follow the periodic variations in the 
mantle's rotation rate. Thus, only the secular variations on the 
mantle will be followed by the core and we may write: coc = 
To — 5, where Hj is the averaged rotation of the mantle and 



6 ^ T s'"'' dt' ^ — 



(28) 



Just before crossing the spin-orbit resonance we have cjq = 
pn + 2Acl>Ik, and thus 

7 CO = - pn ^cJ^ -6 - pn 

= -Au) - — ^ -Au + [p - E(e)] . (29) 

n K,„ K 

When using the present eccentricity of Mercury and the val- 
ues of K and K estimated in previous sections, we compute 
X - 19.5. Substituting this into expression (|27] ) we obtain a 
probability of capture o f 100% in the 3/2 spin- orbit resonance. 
However, as noticed by [Peale and BossI (1977), if we compute 
the probability for the 2/1 resonance we also get 100%. Thus, 
either the planet started its rotation below the 2/1 resonance, 
which is unlikely, or there must be another mechanism to avoid 
capture in the 2/1 and higher order resonances. 

3.3. Effect of planetary perturbations 

The orbital eccentricity of Mercury undergoes important sec- 
ular perturbations from the other planets (FiglH) and its contri- 
bution needs to be taken into account. The mean value of the 
eccentricity is e = 0.198, slightly lower than the present value 
e ^ 0.206, but we also observe a wide range for the eccen- 
tricity variations, from nearly zero to more than 0.45 (Paper I, 
Laskarl |2Q08j). iven if some of these episodes do not last for 
a long time, they will allow additional capture into and escape 
from spin-orbit resonances. Moreover, the capture probabilities 
are also modified for different eccentricities: for the same reso- 
nance we can have zero or 100% of captures depending on the 
eccentricity value (Figl2]). For all resonances, the capture prob- 
ability is 100% whenever the eccentricity is close to the equi- 
librium value for the rotation rate, E(e) - N(e)/Q.(e) (Eq 12111. 
but it tends to decrease as the eccentricity moves away from 
this equilibrium value. If the eccentricity is too high (or too low 
if the spin is increasing from lower rotation rates) some reso- 
nances cannot be reached and the probability of capture sud- 
denly drops to zero (Fig|2]l. 
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Table 2: Critical eccentricity for the resonance p. If e < e^, the resonance 
p becomes unstable, and the solution may escape the resonance (Paper I). The 
critical eccentricity is obtained when jS//(p, e) < K[Q.(e)p - N(e)]. 
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For a non-constant eccentricity e(f), the limit solution of the 
rotation rate when D(j) - (EgfTTli is no longer w/n = E{e), 
but more generally: 



^(0 Jo 
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g{t) = exp 



K r' 

c,„n Jo 



Q(e(T)) t/r 



(30) 



(31) 



The dissipation torques can thus drive the rotation rate sev- 
eral times across the same spin-orbit resonance, increasing the 
chances of capture. 

Another important consequence of a non-constant eccentric- 
ity is that all resonances but the 1/1 may become unstable. In- 
deed, the amplitude of the libration torque depends on the co- 
efficient H{p, e) (Eq|3]l, which goes to zero with the eccentric- 
ity, except for the 1/1 resonance (Tab[Tll. Whenever the am- 
plitude of the libration restoration torque becomes smaller than 
the amplitude of the dissipation torque, equilibrium in the spin- 
orbit resonance can no longer be sustained and the resonance 
is destabilized. Critical eccentricities for each resonance are 
listed in Table |2] obtained when the torques become equivalent 
(Eqsiaa): 



Q.{e)p - Nie) 
H(p,e) 



K 



(32) 



Table 3: Capture probabilities in several spin-orbit resonances (in percentage). 
In the left panel (T only) we consider the effect of tides alone, while in the 
right panel (T + CMF) both tides and core-mantle friction effects are taken into 
account. The first column (Pcap) refers to the theoretical estimation given by 
expression (20), while the next column (num.) refers to the estimation obtained 
running a numerical simulation with 2000 close initial conditions, dift'ering by 
;r/1000 in the libration angle. Planetary perturbations are not considered and 
we used a constant eccentricity e = 0.206. 





T only 


T-t-CMF 
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p 

^ cap 


num. 


p 

-* cap 


num. 




(%) 


(%) 


(%) 


(%) 


5/1 






1.6 


0.3 


9/2 






3.1 


1.3 


4/1 


0.1 




5.9 


4.8 


111 


0.1 


0.1 


11.4 


10.9 


3/1 


0.3 


0.4 


22.6 


22.8 


5/2 


0.7 


1.4 


46.6 


46.2 


2/1 


1.8 


1.7 


100.0 


100.0 


3/2 


7.7 


7.2 


100.0 


100.0 



4.1. Simulations without planetary perturbations 

Before considering the effect of planetary perturbations we 
can test numerically the theoretical estimates of the capture 
probability given by expressions ( |22] ) and ( l27l i. Since cap- 
ture in resonance is a statistical process we need to perform 
many integrations with slightly different initial conditions. For 
that purpose we ran 2000 simulations using a fixed eccentric- 
ity (e - 0.206), initial rotation period of 2 days, zero obliq- 
uity and different initial libration phase angles with step-size 
of ;r/1000 rad. Results are listed in Tabled We can see that 
there is a good agreement between the theoretical previsions 
and the numerical estimation of the probabilities. As discussed 
in sections 13.11 and 13.21 the probability of capture when con- 
sidering only the effect of tides is very small (~ 7% for the 
3/2 resonance), while it becomes very important when core- 
mantle friction is added (100% for the 3/2 resonance). They 
are als o in conformity wit h those obtai ned in the previous stud - 
ies by iGoldreich and Peale C1967.) and iPeale and BossI (Il977h . 
As they all noticed, when the effect from core-mantle friction is 
considered, the probabilities of capture are greatly enhanced for 
all spin-orbit resonances. In particular, capture in the 2/1 res- 
onance also becomes 100%, preventing a subsequent evolution 
to the 3/2 resonance. 



4. Numerical simulations 

We will now use the dynamical equations established in sec- 
tion |2]to simulate the final evolution of Mercury's spin by per- 
forming massive numerical integrations. The main goal is to 
illustrate the effects described in section [3] in particular the 
probabilities of capture and escape from spin-orbit resonances. 
Mercury geophysical models and parameters in use are those 
listed in section |2] We recall here the most uncertain values: 
k2 = 0.4, e = 50 and V = ^-''m^s-i. 



4.2. Inclusion of planetary perturbations 

When planetary perturbations are taken into account, the ec- 
centricity presents chaotic variations with many excursions to 
higher and lower values than today (Laskar, 1990, 1994, 2008, 
Paper 1). It is then impossible to know its exact evolution at the 
time the planet first encountered the spin-orbit resonances. A 
statistical study with many different orbital solutions is the only 
possibility to get a global picture of the past evolution of the 
spin of Mercury. In Paper 1 we performed such a study by in- 
tegrating 1000 orbits over 4 Gyr in the past starting with very 
close initial conditions. This statistical study was only made 
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possible by the use of the averaged equations for the motion of 
the Solar System (Laskar, 1990, 1994). 

In figure [3] we show five examples of the eccentricity evo- 
lution through the 4.0 Gyr We choose some cases illustrative 
of the chaotic behavior, where we can see that the eccentricity 
can be as small as zero, but it can also reach values as high as 
0.5. In some cases the eccentricity can remain within [0.1, 0.3] 
throughout the evoluti on, whi l e in o ther cases it can span the 
whole interval [0, 0.5] (lLaskaiil2008h . 

Owing to the chaotic evolution, the density function of the 
1000 solutions over 4 Gyr is a smooth functio n (Fig|^, wel l 
approximated by a Rice probability distribution (lLaskaiil2008h . 
The eccentricity excursions to higher values allow the planet 
to cross the 3/2 resonance several times, and thus increase the 
probability of capture. This behavior becomes very important if 
the evolution is driven by tidal friction alone. Even though the 
probability of capture in a single crossing of the 3/2 spin-orbit 
resonance is only around 7%, multiple crossings increase it up 
to 55% (Paper I). 

4.3. Planetary perturbations with core-mantle friction 

In presence of an efficient core-mantle friction the multiple 
crossings of the 3/2 resonance are no longer needed, since the 
capture in this resonance after a single crossing is already 100% 
(Fig|2]). Nevertheless, eccentricity excursions to lower values 
can destabilize the equilibrium in any spin-orbit resonance dif- 
ferent from the 1/1 (Tab|2|. This effect was already present 
when the core-mantle friction was not considered (Paper I), but 
with small influence on the results, while here it becomes of 
capital importance. Indeed, it may allow the evasion from pre- 
vious captures in higher order resonances than the 3/2 and per- 
mit subsequent evolution to the present observed spin state. 

In order to check this new scenario, we have performed a sta- 
tistical study of the past evolutions of Mercury's orbit, with the 
integration of the same 1 000 orbits over 4 Gyr in the past used in 
Paper I. We now additio nally include the effect of co re-mantle 
friction as described by iGoldreich and Peald (119671) . i.e., we 
will consider the full dynamics of the spin governed by Eq.(fT3]l. 

Assuming an initial rotation period of Mercury of lOh, we 
estimated that the time needed to de-spin the planet to the slow 
rotations would be about 300 million years. We will then start 
our integrations already in the slow -rotation regime, with a rota- 
tion period of 10 days (lo ^ 8.8«), zero obliquity and a starting 
time of -4 Gyr, although these values are not critical. In TableH] 
we show the amount of captures for each resonance at the end of 
the simulations (column "final"). We also list the resonances in 
which the spin was first captured before being destabilized (col- 
umn "F'cap.") and we recall the results obtained for a constant 
eccentricity (e - 0.206) and in Paper I, with a model without 
core-mantle friction. 

After running 1000 trajectories we observe that the spin of 
Mercury preferably chooses one of the three final configura- 
tions: 5/2, 2/1 or 3/2 (Tabl4]i. With 26% of captures, the present 
configuration no longer represents the most probable final out- 
come, as it was in absence of core-mantle friction (Paper I). 
However, it is still among the most probable scenarios, the al- 
ternatives receiving comparable amounts of captures (22% and 



Table 4: Capture probabilities in several spin-orbit resonances (in percentage). 
We performed a statistical study of the past evolutions of Mercury's spin, with 
the integration of 1000 orbits over 4 Gyr, a initial rotation period of 10 days 
and zero obliquity. In the "l"cap." column we list the resonances in which the 
spin was first captured (before being destabilized). In the "final" column we 
list the results after the full 4 Gyr of simulations. For comparison we also list 
the results obtained with a constant eccentricity e = 0.206 ("const.") and the 
final results obtained in Paper I ("C&L04" ), with a model without core-mantle 
friction. 







number of captures 
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const. 


l^cap. 


final 


C&L04 




(%) 


(%) 


(%) 


(%) 


6/1 




0.1 






11/2 




0.4 






5/1 


0.3 


1.3 






9/2 


1.3 


2.7 






4/1 


4.7 


5.3 






111 


10.3 


8.7 


4.7 




3/1 


19.0 


15.5 


11.6 




5/2 


29.8 


26.5 


22.1 




2/1 


34.6 


31.2 


31.6 


3.6 


3/2 




8.1 


25.9 


55.4 


1/1 




0.2 


3.9 


2.2 


none 






0.2 


38.3 



32% respectively for the 5/2 and the 2/1 resonances). The 5/2 
and the 2/1 spin-orbit resonances benefit from the fact that the 
planet must cross them first. On the other hand, the 3/2 reso- 
nance is more stable and the chances of capture are higher when 
crossed. 

Since the eccentricity of Mercury 4.0 Gyr ago can be around 
0.4 (e.g. Fig|3]l, at the moment of the first encounter with the 
spin-orbit resonances, capture in resonances as high as the 6/1 
can occur (Fig|2|. Because the probability of capture is small 
and because there are not many orbital solutions reaching such 
high values for the eccentricity, we only count about 10% of 
captures in resonances above or equal to the 4/1 (TabHj). Once 
captured, these equilibria can be maintained as long as the ec- 
centricity remains above the respective critical values (Tabl2|i. 

Contrary to the results predicted for a constant eccentricity 
(TablH, we also registered a few trajectories directly captured 
in the 3/2 resonance just after the first passage through the reso- 
nance area. When we used the present value of the eccentricity 
(e = 0.206), the 3/2 resonance could not be attained because 
for that value capture probability in the 2/1 resonance is 100% 
(FiglJl). However, for eccentricity values lower than about 0. 19, 
the probability of capture in this resonance decreases, as well as 
for higher order resonances. For instance, when the eccentricity 
is 0.09, capture in the 2/1 resonance drops to 50%. Thus, since 
the eccentricity of Mercury is varying, it may happen that about 
4.0 Gyr ago its value was much lower than today and the spin 
managed to avoid all the spin-orbit resonances higher than the 
3/2 and was directly captured in the present observed config- 
uration. We estimate nevertheless that the probability for this 
scenario to occur is very low, only about 8% (Tab|4]l. 
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4.4. Critical eccentricities 

Over 4.0 Gyr of evolution the eccentricity has many chances 
of experiencing a period of very small values (e.g. Fig|3]l. Even 
when a period of low eccentricity does not last for a long time, a 
single passage of the eccentricity below a critical value (Tab|2| 
can be enough to destabilize the corresponding spin-orbit reso- 
nance. 

All orbital solutions were generate d starting frorn initial con- 
ditions close to the present values (see lLaskan,l2008l) . and there- 



fore converge to the same final evolution. The eccentricity 
behavior is thus identical for the last 50-60 Myr (FigdJ, be- 
fore which the chaotic diffusion dominates (Fig|5]l. During the 
last 50 Myr the eccentricity certainly reached values lower than 
0.13, thus the 4/1 and above spin-orbit resonances cannot rep- 
resent a possible final outcome for Mercury (64/1 ^ 0.136). For 
the 7/2 and lower order spin-orbit resonances, capture until the 
present day is not forbidden by the last 50 Myr of Mercury's 
evolution, but depends on the true orbital evolution of the ec- 
centricity (Fig|5]l. The higher is the critical eccentricity, the 
lower is the probability of remaining trapped, because more or- 
bital solutions will come below this value. 

In figure|6]we plot the cumulative distribution of the minimal 
eccentricities attained for each one of the 1000 orbital solutions 
that we used. We also mark with straight lines the critical values 
of the eccentricity for each spin-orbit resonance (Tab|2]l and use 
dots to represent the amount of captures obtained numerically 
for spin-orbit resonances that are still stable below each critical 
value of the eccentricity. Since a large amount of the orbital 
solutions experience at least one episode with an eccentricity 
below 0.05 (logjoS ~ -1-3), about 84% of the final evolutions 
will end in the 5/2 spin-orbit resonance or lower (Tab|4]i. By 
comparing the eccentricity instability thresholds for each spin- 
orbit resonance with the amount of captures obtained numeri- 
cally below that resonance we see that there is a good agree- 
ment, suggesting a strong correlation between the orbital evo- 
lution of the eccentricity and the percentage of captures in each 
resonance. The reason why there is not full agreement between 
the two is because spin-orbit resonances below critical values of 
the eccentricity can also be attained by trajectories that escaped 
capture in higher order resonances, that is, they can be attained 
even if the eccentricity is never below the critical value for that 
resonance (Fig|2]i. 

When comparing the results after 4.0 Gyr with those after 
the first capture, we verify that the 5/2 resonance (and above) 
lose a significant amount of previously captured solutions. The 
amount of orbits captured in the 2/1 resonance remains roughly 
the same, because the number of trajectories quitting this reso- 
nance is more or less compensated by the incoming trajectories 
from higher order resonances. The 3/2 is the real winner of this 
transition process, as the amount of trajectories that end in this 
last configuration is about 4 times larger than it was initially. An 
identical scenario was already observed in Paper I, except that 
only a few captures occurred in spin-orbit resonances higher 
than the 2/1 and they were all subsequently destabilized (Tab|4]i. 

As in Paper I, we also noticed about 4% of the trajectories 
captured in the 1/1 spin-orbit resonance (TabHJi. Since the prob- 
ability of capture in the 3/2 spin-orbit resonance is almost 100% 



even for very low values of the eccentricity (Fig|2li, the major 
possibility of evolving into the 1/1 resonance is by destabilizing 
the 3/2. This becomes a possibility if the eccentricity is almost 
zero, that is, for e < 3 x 10 (Tab©. 

4.5. Different scenarios of evolution 

The critical eccentricity needed to destabilize the 2/1 spin- 
orbit resonance is 62/1 ~ 0.0046 (Tab|2]i. Whenever the orbital 
eccentricity is below this value, the spin will then evolve to- 
wards the 3/2 resonance or below (Fig|7]i. However, this is not 
the only possibility of achieving this last configuration if the 
planet was first captured in a higher-order resonance. Indeed, 
as discussed for the V capture column (Tab|4|i, if the eccen- 
tricity is lower than 0.19 at the time the planet crosses the 2/1 
resonance (Fig|2ll, the chances of capture are lower than 100%, 
opening some space for subsequent evolution to the 3/2 reso- 
nance. For instance, for a previous capture in the 5/2 resonance, 
an eccentricity of 65/2 ~ 0.025 will destabilize it and produce 
a capture probability of only 14.4% in the 2/1 resonance, i.e., 
when the 5/2 resonance is destabilized there is about 85% of 
chance of ending in the 3/2 present configuration (Fig|7]i. 

In order to exemplify the multitude of possible evolutionary 
scenarios, we performed another kind of experiment. Adopt- 
ing a particular orbital solution, which presents a gradual de- 
crease in the eccentricity (Figl3^), we integrated close initial 
conditions for the spin. Since for this orbital solution the ec- 
centricity is high at the time of the first encounter, there is a 
great chance of capturing the spin in a spin-orbit resonance with 
p > 5/2. In figure |7] we plot the behavior of four trajectories, 
each one initially captured in a different spin-orbit resonance, 
when the eccentricity approaches a zone of very low values. 
As expected, the spin-orbit resonances are sequentially aban- 
doned as the eccentricity assumes small values. In particular, 
we observe that the resonances are quit immediately after the 
eccentricity is below the critical values listed in Table |2] After 
being destabilized, the spin can evolve directly to the present 
3/2 configuration, or can be trapped in an intermediate spin- 
orbit resonance. In this example, the eccentricity become lower 
than 62/1 ~ 0.0046 around -1.79 Gyr and captures in the 2/1 
resonance become destabilized after that date. 

We purposely plot one situation, where the spin does not end 
in the 3/2 resonance, however In this case, at the moment the 
eccentricity becomes lower than 62/1, the spin is still captured 
in the 5/2 resonance. This resonance logically becomes desta- 
bilized and the rotation rate decreases. Nevertheless, at the mo- 
ment the spin encounters the 2/1 resonance the eccentricity is 
again higher than e2/\, and therefore there is a chance of cap- 
ture in this resonance, preventing a subsequent evolution toward 
the 3/2 state (FigQ. 

4.6. Constraints on the orbital evolution 

We have seen in previous sections that there is an important 
correlation between the minimal eccentricity attained by Mer- 
cury through its orbital evolution and the probability of capture 
in a given resonance (Fig|6]l. The lower is the minimal eccen- 
tricity, the higher is the probability of achieving a low order 
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Table 5: Capture probabilities in spin-orbit resonances (in percentage), when 
the eccentricity descends below a given critical value (Tab[3- 



p 


^3/2 


^2/1 


^5/2 


^3/1 


^7/2 


^4/1 


111 


— 


— 


1.2 


2.6 


3.8 


4.7 


3/1 




2.0 


4.2 


7.2 


10.7 


11.6 


5/2 




3.2 


5.2 


16.0 


20.7 


22.1 


2/1 


3.1 


5.4 


23.9 


30.5 


32.6 


31.6 


3/2 


25.0 


73.2 


55.1 


37.2 


27.7 


25.9 


1/1 


68.8 


15.4 


9.7 


6.1 


4.2 


3.9 


none 


3.1 


0.7 


0.5 


0.3 


0.2 


0.2 



spin-orbit resonance. In particular, each time the eccentricity 
descends below a given critical value for a spin-orbit resonance 
(Tab|2]i, the spin will evolve into a lower resonance. 

Since we know the distribution of the minimal eccentricities 
(Fig|6]l, we can estimate the probability of ending in a specific 
spin-orbit resonance given the value of the minimal eccentric- 
ity of a considered orbit, /"cap/e, ■ For that purpose we eliminate 
all the trajectories for which the minimal eccentricity is above 
the critical value (Tab|2|l and then count the number of captures 
in each resonance for the remaining orbital solutions. Results 
are listed in Table |5] While for an arbitrary orbital solution the 
probability of capture in the present 3/2 spin-orbit resonance is 
only 25.9%, this value rises to 55.1% if we assume that the ec- 
centricity of Mercury was below 65/2 ~ 0.025 at some time in 
the past, or even up to 73.2% if the eccentricity descends below 
62/1 ~ 0.0046. Results for the critical eccentricity e\j\ » 0.136 
are the same as the global results shown in Table |4] because 
the eccentricity of Mercury was below that value in the most 
recent 50 Myr, where the chaotic behavior is not significant 
(FigHl. Notice also that there are always a few captures left 
in resonances above the corresponding critical value. This can 
be explained by the same effect described in the last paragraph 
of section |43] and illustrated in figure |7l 

Inversely, since we know that the rotation of Mercury is 
presently captured in the 3/2 spin-orbit resonance, we can esti- 
mate the probability for the eccentricity to have descended dur- 
ing its past evolution below a specific critical level, Pg^- Using 
conditional probabilities, we have then 



orb 



(33) 



cap 



where /"cap/e, is the probability of ending in a specific spin-orbit 
resonance given the value of the minimal eccentricity (Tab|5]l, 
PoTh the probability for the eccentricity to reach that minimal 
eccentricity (FiglHl and Pcap the global capture probability in 
the specific spin-orbit resonance (TablD). Results for the 3/2 
spin-orbit resonance are given in Table |6] These probabilities 
for the orbital evolution of the eccentricity are the same as if 
we select only the evolutions that finished in the 3/2 spin-orbit 
resonance and then look at the minimal eccentricity distribu- 
tion. From the above analysis we conclude that there is a strong 
probability that the eccentricity of Mercury reached very low 
values; in particular there is about a 77% chance that it de- 



Table 6: Probability for the eccentricity of Mercury to have descended below a 
specific critical level {Pg^ ) given that its rotation is captured in the 3/2 spin-orbit 
resonance today. Results for each critical eccentricity (Tab|2) are obtained 
from Eq.(33) with Pcup = 25.9% (Tab|4). 



p 






65/2 


63/1 




e4/i 




25.0 


73.2 


55.1 


37.2 


27.7 


25.9 


Porb 


3.2 


27.1 


40.3 


64.5 


99.7 


100.0 




3.1 


76.6 


85.7 


92.6 


98.8 


100.0 



scended below e2/i ~ 0.0046 (but only a 3% chance of going 
below 63/2 ~ 0.00003). 

5. Conclusions 

Du e to the increasing evidence of a molt en core inside Mer- 
cury (iNess et al.L 119741 iMargot et al. [ l2007h . viscous friction at 
the core-mantle boundary is expected and its consequences to 
the spin must be taken into account. An important consequence 
is a considerable increase in the probability of capture for all 
spin-orbit resonan ces; in particular, for the 2/1 and the 3/2 it 
can reach 100% CPeale and Bossl Il977h . Since it is believed 
that Mercury's initial rotation was much faster than today, a 
destabilization mechanism is then required to allow the planet 
to escape from the 2/1 and higher order resonances and subse- 
quently evolve to the present observed 3/2 configuration. 

With the consideration of the chaotic evolution of the eccen- 
tricity of Mercury we show that such destabilization mechanism 
exists whenever the eccentricity becomes smaller than a critical 
value for each spin-orbit resonance (TablH. This mechanism 
was already described in Paper I, but becomes of capital impor- 
tance when core-mantle friction is taken into account. There 
are two main possibilities to evolve into the 3/2 configuration: 

• The eccentricity becomes lower than the critical value for 
the 2/1 spin-orbit resonance (e2/i ~ 0.005) and evolves 
into the 3/2. 

• The eccentricity becomes lower than the critical value for 
a higher order resonance than the 2/1, and then crosses 
this resonance with an eccentricity lower than e < 0.19. 
This allows a non zero probability of escaping the 2/1 res- 
onance, and subsequent evolution into the 3/2. 

The other mechanism of capturing in the 3/2 resonance de- 
scribed in Paper I, consisted in a returning to the 3/2 spin-orbit 
resonance after an increase in the eccentricity. This effect is not 
as important when we take into account core-mantle friction, 
since the most part of the trajectories are captured in resonance 
after a single passage. 

After running 1000 orbital solutions, starting from 4Gyr in 
the past until they reached the present date, the spin of the planet 
was captured in a spin-orbit resonance 99.8% of the time. The 
main resonances to be filled and the respective probability were 
(Tabll: 



5/2 



22.1%, 



2/1 



31.6%. 



3/2 



25.9% , 



(34) 
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Although in this case the present configuration no longer repre- 
sents the most probable final outcome, as it was in absence of 
core-mantle friction (Paper I), it is still among the most proba- 
ble scenarios. 

Moreover, if we assume that at some time in the past, the 
eccentricity of Mercury becomes lower than 65/2 ~ 0.025 or 
62/1 ~ 0.005 respectively, the probability of reaching the 3/2 
spin-orbit resonance rises to 55% and 73% respectively (Table 
|5]l. Given that Mercury is presently trapped in the 3/2 config- 
uration, we can also estimate that the eccentricity of Mercury 
has known at least one period of very low eccentricity during 
its past evolution, with about 86% and 77% of chances of being 
below 65/2 ~ 0.025 and 62/1 ~ 0.005 respectively (Table|6ll. 

The probability of capture in the 3/2 resonance can also be 
increased if the orbital eccentricity experiences more periods 
near zero. This can be achieved if we use direct integration of 
the Solar System instead of the averaged equations, because the 
true eccentricity is expected to undergo some additional small 
variations around the value obtained for the averaged equations 
(Laskar, 2008). Alternatively, the probability of capture in the 
3/2 resonance can still be increased if we are able to increase 
the critical eccentricities that destabilize spin-orbit resonances 
(Tab|2]i. This can be achieved if the tidal dissipation is stronger 
(k2 > 0.4 and/or Q < 50) or if C22 < 1.0 x lO^^ (Eql32li. 

Lower values for the core effective viscosity, v, will not 
change the critical eccentricities, but will decrease the amount 
of captures for all spin-orbit resonances. As a consequence, it 
becomes easier to escape from the capture in spin-orbit reso- 
nances, and all those trajectories that also escape the 3/2 reso- 
nance can be later trapped there when the eccentricity experi- 
ences a period with e > 0.325 (Paper I). We then believe that 
the true scenario for the evolution of the spin of Mercury may 
be somewhere between the scenario described here, with an ef- 
ficient core-mantle friction effect, and the scenario described in 
Paper I, for a total absence of core-mantle friction. In the future, 
different dissipative parameters and models could be tested as 
well as the effect of the obliquity, that was supposed to be zero 
in the present study. 
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Figure 1 : Evolutio n of Mercury's eccentricity over 50 Myr in tlie past 
iLaskar et alll2004 jihl) . 



Yoder, C. F., 1997. Venusian Spin Dynamics. In: Venus II: Geology, Geo- 
physics, Atmosphere, and Solar Wind Environment, pp. 1087-1124. 
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Figure 3: Some examples of the possible variations of the eccentricity of Mer- 
cury thi'ough the past 4.0 Gyr The eccentricity can be as small as zero, but 
it can also reach values as high as 0.5. In some cases the eccentricity can re- 
main within [0.1, 0.3] throughout the evolution, while in other cases it can span 
the whole interval [0,0.5]. All these solutions converge to the known recent 
evolution of the planet's orbit (Fig|5). 
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Figure 4: Probability density function of Mercury's eccentricity (Paper 1). Val- 
ues ar e computed over 4 Gyr for the numerical integration of the secular equa- 
tions jLaskar et all l2004a b : Laskar. 2008) for 1000 close initial conditions 
(LA04). The mean value of the eccentricity is e = 0.198. 
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Figure 2: Probability of capture in some spin-orbit resonances for different values of the eccentricity, under the effect of tides and core-mantle friction (Eq |27r The 
dashed line corresponds to a planet increasing its spin from slower rotation rates, while the solid line corresponds to a planet de-spinning from faster rotation rates. 
For all resonances, capture probability is 100% whenever the eccentricity is close to the equilibrium value for the rotation rate, w/n = E(e) (Eq |2U . It suddenly 
decays to zero when the equilibrium rotation rate falls outside the resonance width, i.e., the tidal evolution prevents the planet from crossing the resonance. 
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Figure 7: Four possible final evolutions for the spin of Mercury. Adopting a particular orbital solution, which presents a gradual decrease in the eccentricity (Fig[3ji), 
we integrate close initial conditions for the spin. We observe that spin-orbit resonances are quit immediately after the eccentricity is below the critical values listed 
in Table 12] After being destabilized, the spin can evolve directly to the present 3/2 configuration, or can be trapped in an intermediate spin-orbit resonance. We 
puiposely left one situation where the spin does not end in the 3/2 resonance. At the moment the eccentricity becomes lower than £2/1 ~ 0.0046, the spin is still 
captured in the 5/2 resonance. This resonance is then destabilized, but when the spin encounters the 2/1 resonance the eccentricity is already higher than £2/1 ■ Thus, 
there is a chance of capture in this resonance, preventing a subsequent evolution toward the 3/2 state. 
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Figure 5: Some examples of the recent evolution of the eccentricity of Mer- 
cury. During the last 50-60 Myr all the orbits present the same evolution, before 
which the chaos effect takes place. Horizontal lines correspond to the critical 
eccentricities that destabilize the equilibrium in a given spin-orbit resonance 
(Tab|2). During the last 50 Myr the eccentricity was certainly below 0.13, thus 
the 4/1 resonance is not a possible final outcome for Mercury. For some orbits 
the eccentricity is below 0.09 in the last 100 Myr and the 7/2 resonance was also 
destabilized. However, since this scenario is not true for all orbital solutions, 
we may expect a few final evolutions captured in this last configuration (Tab|4). 
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Figure 6: Cumulative distribution of the minimal eccentricities attained for the 
1000 orbital solutions that we used. Straight lines represent the critical values 
of the eccentricity for each spin-orbit resonance (Tab|2), while dots represent 
the amount of captures obtained numerically for spin-orbit resonances that are 
still stable below each critical value of the eccentricity (Tab|3). 



14 



